An ab-initio calculation of magnetic resonant x-ray scattering spectra in NiO 
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We investigate the magnetic resonant x-ray scattering spectra around the K edge of Ni in an- 
tiferromagnetic NiO using an ab-mitio band-structure calculation based on the density-functional 
theory. By taking account of orbital polarization through the spin-orbit interaction, we reproduce 
well the spectra obtained experimentally, thus demonstrating the usefulness of the ab-initio cal- 
culation. It is shown that the main-edge peak, which mainly comes from the dipolar (Is — > 4p) 
transition, is a direct reflection of the orbital polarization of the 4p states. It is clarified that the 
4p orbital polarization is mainly induced from the spin polarization on the 4p states by the spin- 
orbit interaction. The 3d orbital polarization at neighboring Ni sites gives rise to only a minor 
contribution to the 4p orbital polarization through a p-d mixing. It is also shown that the pre-edge 
peak, which mainly originates from the quadrupolar (Is — > 3d) transition, is a direct reflection of 
the orbital polarization of the unoccupied 3d states. It shows a Fano-type antiresonant dip due to 
interference with the nonresonant contribution, in agreement with the experimental result. 

PACS numbers; 78.70.Ck, 71.20.Be, 75.50.Ee 



I. INTRODUCTION 

Resonant x-ray scattering (RXS) is a powerful tool 
for investigating the properties of magnetic and orbital 
orders. Resonant enhancements on antiferromagnetic 
(AF) superlattice spots have been observed in several 
transition-metal compounds such as RbMnFa (Ref. Q), 
KCuFg (Ref.H, CoO (Ref.lS) and NiO (Refs.ii). RXS 
is described as a second-order process wherein a photon 
is virtually absorbed upon exciting a core electron to un- 
occupied states, and then emitted upon recombining the 
excited electron with the core hole. For the id transition- 
metal if-edge, the intermediate states involve unoccupied 
4p states in the dipolar (El) transition and unoccupied 
?>d states in the quadrupolar (E2) transition. The RXS 
intensities due to the El and E2 transitions constitute 
main-edge and pre-edge structures, respectively. The 
main-edge intensity is sensitive to the electronic struc- 
ture of neighboring sites since the unoccupied 4p states 
are highly extended in space. This characteristic is now 
well recognized in RXS on orbital superlattice spots in 
the K edge of transition-metal compounds. At an early 
stage, a mechanism was proposed, according to which 
the intensity is related to the 4p states modified by the 
intra-atomic Coulomb interaction with the Zd electrons 
constituting the orbital order^^, but extensive studies 
based on band-structure calculation o^i^d" revealed that 
the intensity is related to the 4p states modified by the 
lattice distortion. 

In our previous work,^^ we studied the mechanism of 
magnetic RXS in the antiferromagnetic phase of NiO, 
and analyzed the experiment performed by Neubeck et 
al^ To deal with the extended nature of the 4p states, 
we used a tight-binding model. Since orbitals must cou- 
ple to the magnetic order to generate the intensity of the 
magnetic RXS, it is essential to take account of the spin- 
orbit interaction (SOI). The spin polarization alone in 



the Ap states cannot give rise to the intensity, because, if 
we add the amplitudes due to up-spin electrons to that 
due to down-spin electrons, the total becomes the same 
at each site, leading to a cancellation of amplitudes be- 
tween different magnetic sublattices. Through the model 
calculation, we clarified that the main-edge intensity is 
generated from the orbital polarization in the 4p states. 
Generally speaking, such polarization is brought about 
by the ip SOI, by the intra-atomic Coulomb interaction 
between the Ap and 3d states, and by the mixing of the 
Ap states with the 3d states of neighboring Ni atoms. Al- 
though we obtained the magnetic RXS spectra that were 
in qualitative agreement with the experiment, it is still 
unclear how the Ap states are orbitally polarized. In this 
paper, to eliminate the uncertainty of the model calcula- 
tion, and to clarify the mechanism of the magnetic RXS, 
we carry out an ab-initio calculation of the RXS spectra 
using a band-structure calculation based on the density- 
functional theory. We take account of the SOI to deal 
with the orbital polarization. We neglect the attractive 
interaction between the photoexcited electron and the Is 
core hole. This effect on the Ap states is known to be 
small, from the analysis of the magnetic circular dichro- 
ism (MCD)i2iiii4 and of the RXS on orbital superlattice 
spots at the K edge of transition nietals.^*i2ii^ 

In this paper, we show that the ab-initio calculation 
reproduces well the experimental spectra. Moreover, we 
calculate the spectra by turning the SOI on and off se- 
lectively in several states in order to clarify how the Ap 
orbital polarization is brought about. The result is that 
the main-edge intensity is slightly reduced when the SOI 
on the 3d states is turned off, while it is largely dimin- 
ished when the SOI on the Ap states is turned off. This 
indicates that the 3d orbital polarization has little influ- 
ence on the Ap orbital polarization. The same conclusion 
has been drawn on the basis of an ab-initio calculation 
of the magnetic RXS spectra for KCuFa (Ref. ITsI) and a 
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recent analysis of MCD in MngGaC (Ref.T^. 

We also calculate the contribution of the E2 and non- 
resonant processes. We obtain a single peak in the pre- 
edge region. It is shown to originate from the orbital 
polarization in the 3d states. Although the 3d orbital 
polarization is much larger than the 4p orbital polariza- 
tion, the obtained pre-edge peak is rather small compared 
to the main-edge peak. The peak position is at a higher 
energy than the experimental one. We can ascribe this 
discrepancy to the neglect of the interaction between the 
Is core hole and the 3d electrons in the present calcula- 
tion. We also reproduce well the Fano-type antiresonant 
dip due to the interference between the E2 and the non- 
resonant contributions. 

This paper is organized as follows. In Sec. ^ we show 
the results of band-structure calculations and discuss the 
ground state in the antiferromagnetic phase of NiO. In 
Section llTTI we briefly describe the formulae used in calcu- 
lating the magnetic RXS spectra. In Sec. IIVI we present 
the calculated magnetic RXS spectra and compare them 
with experimental data. Sectionals devoted to conclud- 
ing remarks. 



II. ELECTRONIC STRUCTURE OF NICKEL 
OXIDE 

NiO has the NaCl-type structure with the lattice con- 
stant of a = 4.177 A (Ref. 17). Nickel atoms form an 
fee lattice, as illustrated in Fig. ^ Type-II AF order de- 
velops below Tn = 523 K. The order parameter is char- 
acterized by a wave vector Q, which is directed to one 
of four body-diagonals in the fee lattice. We selectively 
assume that Q = (5^5) in the present study. The mag- 
netic moment lies in the plane perpendicular to the AF 
modulation direction; three directions. Si = [—1,-1,2], 
S2 = [2,-1,-1], or S3 = [-1,2,-1] can be taken for 
Q = (i^^)- Note that the experiment of the magnetic 
RXS in NiO has been performed under the condition that 
Q is specified but domains are formed with respect to the 
spin direction^iii. 

The electronic band structure of NiO is calculated us- 
ing the full-potential linearized augmented-plane-wave 
(FLAPW) methodiSiia in the local-density approxima- 
tion (LDA) . The local exchange-correlation functional of 
Vosko, Wilk and Nusair^° is employed. The experimental 
lattice constant is used and the mufhn-tin (MT) radii are 
set to be 2.2 a.u. for nickel and 1.6 a.u. for oxygen. We 
also perform LDA+U calculation ^iiS^ with the effective 
Coulomb parameter UcS = U — J ior the Ni 3d states, 
in order to investigate the effect of the band gap and the 
magnetic moment on the RXS spectra. 

Table HI shows the calculated values of the band gap 
and the magnetic moments in comparison with the ex- 
periments. LDA calculation considerably underestimates 
the magnetic moment and band gap magnitudes. Intro- 
ducing UcS enhances them but the obtained magnitudes 
are still smaller than those in the experiments. 
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FIG. 1: Sketch of a NiO crystal (only Ni atoms are shown). 
Wave vector Q = (55^) characterizes an antiferromagnetic 
modulation direction. Three arrows denoted as Si, S2, S3, in 
a plane perpendicular to Q indicate possible spin directions 
{S domains). Spins for solid, gray, and open circles align 
alternately. 



TABLE I; Calculated band gap and magnetic moments, com- 
pared with those in the experiments. 25 and L are spin and 
orbital moments, respectively. In the LDA+U calculations, 
the effective parameters of (1) Ucft = 2 eV and (2) Uca = 5 
eV are employed. 





Sgap [eV] 


25' [h] 


L [h] 


L + 2S [h] 


L/S 


LDA 


0.408 


1.210 


0.151 


1.361 


0.25 


LDA+f/(l) 


1.438 


1.467 


0.187 


1.654 


0.25 


LDA+f/(2) 


2.839 


1.647 


0.233 


1.880 


0.28 


Exp." 


4.3 ' 


1.88 


0.32 


2.20, 1.90" 


0.34 



"Unless noted, Ref. I2I 
''Reference |24t 
"^Reference l25t 



Figure 121 shows the density of states (DOS) projected 
on the p and d symmetric states inside the MT sphere at 
a Ni site, calculated in the LDA. The origin of the energy 
is set at the top of the valence band. The right part of the 
figure is shown on a 20-times magnified scale. The p and d 
symmetric states are sometimes referred to as the 4p and 
3d states, respectively, because the wave functions are 
similar to atomic Ap and 3d wave functions inside the MT 
sphere. The magnitude of the 4p DOS is much smaller 
than that of the 3d DOS, because most of the ip density is 
distributed over the interstitial region or on neighboring 
sites. Introducing large t/cff, such as Ucff = 5 eV, in the 
LDA+t/ calculation considerably changes the 3d DOS 
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(not shown), while the Ap DOS remains almost the same. 
Also, exchange splitting in the 4p DOS is negligibly small 
in contrast with the large exchange splitting in the 3d 
states. Both observations suggest that the effect of the 
3d states on the 4p states is very small. 
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FIG. 2: Projected density of states of the Ni 3d and 4p states 
calculated in the LDA. All energies are given with respect to 
the top of the valence band. 



where k, = ki/|k,|, k/ = k//|k/|, and G = G/|G|. 
The quantities S(G) and L(G) in Eq.(|2Jl are the Fourier 
transforms of the spin and orbital densities, respec- 
tivelji^i. Their approximate forms for calculation are 
given in the next section. 
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FIG. 3: Geometry of x-ray scattering. Incident photon with 
wave vector ki and polarization cr or tt is scattered into the 
state with wave vector k/ and polarization a' or tt' at Bragg 
angle Q. The sample crystal is rotated by azimuthal angle t/i 
around scattering vector G = k/ — k^. 



III. FORMULATION FOR MAGNETIC RXS 

We consider the situation that a photon with the en- 
ergy fiw, the momentum k^, and the polarization fj, {= a 
or vr) with the polarization vector e, is scattered into 
the state with the same energy hio (elastic scattering), 
the momentum ky, and the polarization /z' (= a' or tt') 
with the polarization vector e . The scattering geometry 
is depicted in Fig. |31 The magnetic scattering ampli- 
tude is derived from Fermi's golden rule to the second 
orderi2Si2L2Si22i22i^, in which the amplitude consists of 
resonant and nonresonant scatterings. We consider the 
resonant El and E2 transitions and the nonresonant con- 
tribution. The cross section for the elastic scattering at 
a magnetic superlattice spot is expressed as 

^ (X \T^^^> (G, c^) + J^^^, (G, c^) + L^^^, {G,Lu)\\ 

(1) 

where G = k/ — k; is a scattering vector. 

The nonresonant term r^^^'(G,w) is given by 

T,^^'{G, Lo) = (\UG) ■ A" + S(G) . b) , (2) 

where m is the electron mass and c is the light velocity. 
The vectors A" and B are given by 

A" A' - (A' • G)G, A' -4sin2 0{e' x e), (3) 
B = e' X e + (k/ X e')(k/ • e) - (k^ x e){k, ■ e) 

~(k/ X e') X (k, X e), (4) 



The resonant contribution from the El transition is 
given by 

J^^^,{G,uj) = Y,iP"'')c.M^c.'(G,uj)P^,, (5) 

with 

A'W(G,tj) 

~ VN huj-oJAg + ir 

n,A ^ 

(6) 

where = {Ej^—Eg)/h. This represents a second-order 
process wherein a photon is virtually absorbed upon ex- 
citing the Is electron to the 4p states and then emit- 
ted upon recombining the excited electron with the Is 
core hole. The geometrical factors pM p'^ Q^e explic- 
itly written in the Appendix of Ref. [T^. Eg represents 
the ground-state energy. The intermediate state |A) con- 
sists of an excited electron on the 4p states and a hole 
on the Is state with energy E\. N represents the num- 
ber of Ni sites, and n is the index for the Ni sites. The 
dipole operators Xa{n) at site n are defined as xi{n) — x, 
X2(j>-) — y, X3{n) = z in the coordinate frame fixed to the 
crystal axes with the origin located at the center of site n. 
Note that the amplitudes on magnetic spot G — (-Hf) 
are nearly canceled out between different magnetic sub- 
lattices. Only small values remain in an antisymmetric 
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form, 



M(G,w) = 




(7) 



with a, 6, and c being some complex numbers. 

The resonant contribution from the E2 transition is 
given by 



77' 



(8) 



with 



-iGR„ 



n,A' 



12 Vc/ 



\ 2 mt^i, (5 1 (n) I A') ( A' I (n) I (?) 



(9) 



where uja'q = {E\i — Eg)/h. The geometrical factors 

Q^, Q'^ are given in the Appendix of Ref. The 
intermediate states |A') consist of an excited electron 
on the Zd states and a hole on the Is state with en- 
ergy i?A'- The quadrupole operators z^{n) are defined 
as z,{n) = (^/3/2)(a;2 - tf), z^in) = (l/2)(3z2 _ r^), 
2:3(71) = -s/Syz, Z4(n) = \/?,zx, z^{'n) = \^xy. Similar 
to the El transition, the amplitudes on magnetic spot 
G = (ff |) are nearly canceled out between different 
magnetic sublattices, resulting in an antisymmetric form, 



coincides with the experimental one {fku = 8347 eV for 
the peakS.). We set the core- hole lifetime width F as 1 
eV. Figure 0] shows the calculated results and those of 
the if-edge fluorescence experiment.^ The experimental 
spectral shape is reproduced well, suggesting that the 
core-hole potential plays only a minor role in the spec- 
tral shape. As already mentioned in Sec. I, this obser- 
vation is consistent with the recent analysis of the K- 
edge MCD spectra of MngGaC (Ref. [H) as well as of 
if-edge RXS on orbital superlattice spots in transition- 
metal compoundsiSiiS 
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FIG. 4: Absorption coefficient A{uj) around the K edge and 
the experimental fluorescence spectrum (Ref. |^. The core- 
liole energy was adjusted such that the peak position coin- 
cided with that in the experiment. 
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(10) 



with d ^ u being some complex numbers. 

IV. MAGNETIC RXS SPECTRA 

A. Absorption coefficient 

The absorption coefficient within the El transition is 
given by 



A(a;)«^|(AK(j)|5)|^ 



A,a 



(11) 



We neglect the core-hole potential working on the 4p 
states in the final state of the absorption process, and re- 
place I A) by the unoccupied 4p states given by the band- 
structure calculation. The transition matrix elements are 
evaluated using the LDA wave-functions. The core-hole 
energy is adjusted such that the calculated peak-position 



B. Magnetic RXS spectra 

We evaluate the nonresonant term using the "spheri- 
cal approximation" and "dipolar approximation" for the 
Fourier transforms S(G) and L(G), respectively"^^. S(G) 
is expressed as 



with 



S(G) 



(jOr, 



1 



(jOn)(Sr 



Jo{Gr)[<jy^i{r)Yr'dr, 



(12) 



(13) 



where jo(a^) = sinx/x. Rmt is the MT radius and 4>ni{T) 
is the radial function with I being the angular momentum 
on the Ni n site. (S„) is the average of the spin operator 
within the MT sphere at site n. The 3d states (/ = 2) 
give a dominant contribution. 
L(G) is expressed as 



L(G) 



-jG R„ 



(30n)(Ln) 



(14) 
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with 



•J On I 



-Rmt 



(15) 



where goix) = — cosx). (L„) is the average of the 
angular momentum with respect to the angular part of 
the wave function within the MT sphere at the Ni n 
site. Again, the Sd states give the dominant contribution. 
Figure shows calculated |S(G)|///V and |L(G)|/ViV 
and those in the experiment. The calculated values are 
slightly smaller than the experimental ones. This is be- 
cause the LDA calculation underestimates the magnetic 
moments. The Fourier transform S(G) can be directly 
evaluated without the spherical approximation from the 
LDA wavefunctions. We compare the results with those 
evaluated from the spherical approximation. The values 
in units of h are |S(G)|/ViV = 0.57(0.56), 0.32(0.33), 
and 0.10(0.13) for G = {\\\), (|||), and (||f), re- 
spectively, where the numbers in parentheses are values 
obtained by the spherical approximation. The spherical 
approximation works rather well for small values of |G|. 
A similar tendency is expected to exist for L(G) calcu- 
lated with the dipolar approximation, according to the 
study for Ni++ by Blumc.-^i 
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FIG. 5: Form factors of the spin and orbital moments. Cir- 
cles: the spin moment per Ni site, |S'(G)| = \S{G)\/\fN. 
Squares: the orbital moment per Ni site, |L'(G)| = 
|L(G)|/ViV. 

Once S(G) and L(G) are obtained, the nonresonant 
contribution is easily evaluated. Then the resonant con- 
tribution is evaluated by replacing the intermediate state 
|A) in Eq. (3.8) and |A') in Eq. (3.11) with the unoccupied 
p and d symmetric states given by the band-structure 
calculation. We confine ourselves to the calculation on 
a single domain of Si with the AF modulation vector 
Q fixed at (5^5). FigureEIa) shows the spectra in the 
a ^ n' channel for G = (5I5) a function of photon 
energy. The lower panel of Fig. |H{a) shows the calcu- 
lated spectra at the azimuthal angle ip = 270°, where ip 
is defined such that the scattering plane at ■0 = con- 
tains the [1,-1,0] crystal axis. The core- hole energy has 



already been determined such that the calculated peak 
in the K-edge absorption coefficient coincides with the 
experimental one. The upper panel of Fig. El^a) shows 
the experimental curve at ^ = 0°, where the absorption 
correction is made.^''^ The experiment is carried out on 
a sample in which three domains of Si, S2, and S3 are 
mixed with various relative volumes. 

We obtain several peaks at around hoj — 8347 eV 
due to the El transition. The calculated curve agrees 
well with the experimental one. Such comparison makes 
sense, because the spectral shape is a function of pho- 
ton energy and depends little on domains and azimuthal 
angle. The El intensity arises from the orbital polar- 
ization in the 4p states. To clarify how the polarization is 
induced, we calculate the spectra by turning the SOI on 
and off separately on the 3d states and on the 4p states. 
The results are shown in the upper panel of Fig. IHI^b). 
The former slightly suppresses the El intensity, while 
the latter does so strongly. This indicates that, for the 
4p orbital polarization, the contribution of the 3d orbital 
polarization through the p-d hybridization is quite small 
and the 4p orbital polarization is mainly induced by the 
Ap spin polarization through the SOI. This is consistent 
with the recent analysis of the K-cdge MCD in the ferro- 
magnetic phase of MnsGaC (Ref. 16) where it was clari- 
fied that the Ap orbital polarization near the Fermi level 
is mainly induced by coupling to the 3d orbital polariza- 
tion, while that away from the Fermi level is induced by 
the Ap spin polarization through the SOI. 

In addition to the main-edge peak, we obtain another 
peak at around hu = 8333 eV. This pre-edge intensity 
is attributed to the E2 transition, that is, the transition 
to unoccupied 3d states in the intermediate state. As 
shown in the upper panel of Fig. IHl^b), partially turning 
the SOI off generates small El intensities around hut = 
8333 eV, but the net contribution from the El transition 
(the solid line in the upper panel of Fig. El^b)) becomes 
negligible due to a cancellation. As shown in the lower 
panel of Fig.EJb), the E2 intensity is strongly suppressed 
by turning the SOI off on the 3d states, indicating that 
it originates from the 3d orbital polarization. The mag- 
nitude of the pre-edge peak is rather small compared to 
that of the main-edge peak, though the 3d orbital polar- 
ization is much larger than the 4p orbital polarization. 
The experimental pre-edge peak is located at huj — 8330 
eV, which is 3 eV lower than the calculated position. 
We can attribute this difference to the effect of the in- 
teraction between the core hole and the 3d states. The 
attractive interaction is expected to be stronger on the 
3d states than on the Ap states, because the former states 
are more localized than the latter ones. The nonresonant 
contribution is extremely small at this azimuthal angle. 

In Fig. [7| we show the RXS spectra calculated with 
the LDA and the LDA+U. The intensity is enhanced 
due to the local Coulomb interaction UcS- This is con- 
sistent with the enhancement of the orbital moment (see 
Table I). The pre-edge peak is shifted to higher energy 
on introducing UcS, which correlates with the opening 
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larger energy gap. Nevertheless, the effect of UcS on the 
whole spectral shape is small. 

Figure |H1 shows the spectra in the a ^ a' channel 
for G = (55^). Both the calculated and experimen- 
tal curves are for ip = 0. The absorption correction is 
not made on the experimental curve, so the nonresonant 
intensities below and above the resonant peak are not the 
same in the experimental curve^. Since the El transition 
is forbidden for all tp in the a a' channel, no main-edge 
peak appears. The pre-edge peak is located at an ~ 3 
eV higher position, due to the same reason as in the case 
of the a ~* tt' channel. We obtain a Fano-type antireso- 
nant dip on the low-energy side of the peak, which arises 
from an interference with the nonresonant contribution, 
in agreement with the experiment. 

Finally we calculate the azimuthal-angle dependence 
of the intensity. Although the calculation is performed 
for a single domain and it cannot be directly compared 
with the experimental results, it may serve as a guide for 
analyzing mixed domains. Figure El shows the result for 
G = (555)- The upper (lower) panel shows the intensity 
of the main-edge (pre-edge) peak at huj = 8347 eV (8333 
eV) . The dependence of the main-edge peak is controlled 
by the geometrical factor of the El transition. As already 
mentioned, the El transition is completely forbidden in 
the (T — > (t' channel. 



V. CONCLUDING REMARKS 

We presented a theoretical study of the magnetic RXS 
around the K edge of Ni in antiferromagnetic NiO using 
the FLAPW band-structure calculation in the density- 
functional theory. The calculation reproduced the ex- 
perimental spectra, thus demonstrating the usefulness of 
the ab-initio calculation. The main-edge peak originates 
from the orbital polarization in the Ap states through the 
El transition. We found that the 4p orbital polarization 
is mainly generated from the 4p spin polarization through 
the SOI and the effect of coupling to the 3d orbital polar- 
ization at neighboring Ni sites is found to be quite small. 



The 4p orbital polarization is closely related to the K- 
edge MOD spectra. It is known from the analysis of the 
MOD spectra in MngGaC (Ref.|l3), and in Fe, Co and Ni 
(Ref. that the MCD signals near the Fermi level are 
generated mainly by the 4p orbital polarization coupled 
to the 3d orbital polarization at neighboring sites, while 
the MCD signals well above the Fermi level are generated 
from the 4p orbital polarization coupled to the 4p spin 
polarization through the SOI. The present finding for the 
main-edge peak is corresponding to the latter. 

We also obtained the pre-edge intensity through the E2 
transition. The contribution of the El transition to the 
pre-edge intensity is found to be negligible. This shows a 
contrast with the RXS spectra as to the orbital-ordering 
spots in LaMnOs. The El transition generates pre-edge 
intensity through mixing 4p states with neighboring 3d 
states. In LaMnOs the pre-edge intensity of the El tran- 
sition is much smaller than the main-edge intensity but 
still larger than the pre-edge intensity of the E2 transi- 
tion. In the present case, the main-edge intensity on the 
magnetic spot is much smaller than that on the orbital- 
ordering spot in LaMnOa, and so becomes the pre-edge 
intensity of the El transition. Thereby the intensity of 
the E2 transition becomes relatively large. 

The calculated pre-edge peak position relative to the 
main-edge peak is ~ 3 eV higher than the experimental 
one. This suggests that the core-hole potential is more 
effective on the 3d states. The inclusion of such an ef- 
fect by extending the calculation scheme is left for future 
studies. 
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FIG. 6: (color online) (a) Magnetic scattering spectra in the 
o" — > tt' channel for G = (|^^), as a function of photon en- 
ergy. The upper panel shows the experimental spectrum at 
'0 = 0° with three domains mixed and the absorption cor- 
rection is performed^ The lower panel shows the calculated 
spectrum nt ip — 270° for the domain of Si. Resonant El, 
E2, and nonresonant contributions are included in this calcu- 
lation, (b) Magnetic scattering spectra by turning the SOI 
off separately on the 3d states and the 4p states, respectively. 
Nonresonant contribution is not included in these spectra. 
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FIG. 7: Comparison between the RXS spectra calculated 
by the LDA (sohd line) and the LDA+!7 with Ues = 2 eV 
(dashed line). 
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FIG. 8: Magnetic scattering spectra in the ct — > cr' channel for 
G = (iil),asa function of photon energy. The upper panel 
shows the experimental spectra at i/) = without absorption 
correctioni- The lower panel shows the calculated results at 

i' = o. 




FIG. 9: Azimuthal-angle dependence of the magnetic scat- 
tering intensity for G = (Hf)- The upper panel is for the 
pre-edge peak at 8333 cV (non-resonant contribution is ehm- 
inated). The lower panel is for the main peak at 8347 eV. 



